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Abstract 



This second part deals with applications of a general method to describe the quan- 
tum time evolution determined by a Schrodinger equation with time-dependent Hamil- 
tonian. A new aspect of our approach is that we find all solutions starting from one 
special solution. The two main applications are reviewed, namely the Bloch equations 
and the harmonic oscillator with time-dependent frequency. Even in these well-known 
examples some new results are obtained. 

in 
o . 

§ ; 1 Introduction 

The study of quantum models with a time-dependent Hamiltonian is relevant for many 
branches of physics. Instead of solving the Schrodinger equation for a quantum system 
together with its environment, it is often possible to replace the action of the environment 
by a time- dependent term in the Hamiltonian of the system. A slightly more general problem 
is that of solving the von Neumann equation 

^p(t)=i\p(t),H(t)), (1) 
which describes the time evolution of the density matrix p. Its formal solution is 

P (t) = u(t) P (o)u(ty. (2) 

The unitary operators U (t) satisfy the equation of motion 

i(^U(t))u*(t) = H(t). (3) 
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A considerable body of knowledge exists about solving these equations in the case that 
both the Hamiltonian and the density matrix (minus an operator commuting with the Hamil- 
tonian) are linear combinations of the generators Si, S2, ■ ■ ■ , S n of a finite Lie algebra. This 
knowledge has been reviewed in the first part of the present work [31] , hereafter called Part I. 
In addition, a method was presented to obtain expressions for U (t) in a systematic manner. 
The aim of what follows is to show that the method is capable of reproducing known results 
and of obtaining new results. 

Sections 2 and 3 of the paper deal with the SU(2) symmetry. Section 2 discusses the 
Bloch equations and a generalisation involving Jacobi's elliptic functions. Section 3 discusses 
the case when phase modulation is included. Section 4 deals with SU(1,1) symmetry, more 
specifically, the harmonic oscillator with time-dependent frequency. Section 5 considers more 
general oscillators. Finally, Section 6 contains a short discussion of the obtained results and 
of the possibilities for further work. 

2 SU(2) 

The generators of SU(2) satisfy the commutation relations 

[Si, £2] = iSs and cyclic permutations. (4) 

A well-known application of this Lie algebra concerns the Bloch equations — see [3]. It is 
treated in the present Section. The modification obtained by considering phase modulation 
is treated in the next Section. The combination of amplitude and of phase modulation was 
considered in [5] but will not be considered here. 

Note that the incorporation of the Bloch equations into the Maxwell-Bloch equations 
have been studied by many authors, including [TDl [26] . 

More general applications of SU(2) symmetry have been considered in the literature as 
well. Campolieti and Sanctuary [5] applied the Wei-Norman technique to field modulation 
in NMR. Zhou and Ye [12J study the case where all coefficients are time-dependent. They 
introduce Euler angles with the same purpose as in the present work. Finally, Dasgupta [17] 
studies the Jaynes-Cummings model with time-dependent coupling between the spin and the 
photon field. 

2.1 The Bloch equations 

A magnetic spin in a magnetic field is usually described by a Hamiltonian of the form 

H = -ea 3 - lf<7i. (5) 

The Pauli matrices 0-1,0-2,03 are related to the generators of the Lie algebra by o a = 2S a . 
One has h = (— £, 0, e) T . The equations of motion a = hxa are known as the Bloch equations 
— see [I]. Written explicitly, they are 

di = — ea 2 , 
d 2 = eai + £a 3 

d 3 = -£a 2 . (6) 
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2.2 Special solutions 

Assume that the constant e does not depend on time. If £(£) does not depend on time then 
any solution x(t) of the harmonic oscillator equation x + (e 2 + £ 2 )x = determines a solution 
°f ©, given by 

o - ccr, (7) 





ex 


ai = 




a 2 = 








a 3 = 


e 2 + e : 



with integration constant C . 

When is time-dependent then a solution is known only in very specific cases. One 
such case involves Jacobi's elliptic functions sn, cn, and dn, with elliptic modulus k. Let 

= 2ujkcn(ut-k). (8) 

Then a solution of ([6]) exists of the form 

a,i(t) = e cn(ut; k) , 

a,2{t) = ijJ sn(ut; k) dn(u;t; k), 

^{t) = — coksn 2 (ujt; k) + 7, (9) 

with 

f J 2 

7 = - L s = ^-- (10) 

In the limit k = 1 these expressions lead to the well-known result (see Eq. 4.21 of [I]) 

£(t) = 2cj sech(u;t), 
ai(t) = esech(a;t), 
a 2 (t) = co>tanh(u;t) sech(a;i), 

a 3 (t) = -cjtanh 2 (cut) + 7. (11) 
When cl> = e then the limit = can be taken. The rather trivial result is 

m = 0, 

ai(t) = ujcos(idt), 
ct 2 (t) = usm(u;t). 

a 3 (t) = 0. (12) 



2.3 Automorphisms 

Consider the special solution ([9]). Following the general method of Section 4.4 of Part I the 
transformation V(t) is determined by two angles <j>(t), 6{t). The special solution a(t) at time 
t = reads 

a(0) = (e,0, 7 ) T . (13) 
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It is rotated into the fixed vector A(0, 1, 0) T . One has 

z(t) = \J a\(t) + a 2 2 {t) = sjt 2 cn 2 (uot] k) + uj 2 sn 2 (cut; k) dn 2 (wt; k), 

z{0) = e 

A = y/z*(t) + 4(t) = ^+7- (14) 

The angles <p(t) and 9(t) are determined by (62) of Part I. In particular, at t — is 

Bin(0(O)) = ^ = 1 and cos(0(O)) = = 0, 

z(0) z(0) 

sin(0(O)) = -^4^ = -? and cos(0(O)) = ^ = (15) 

A A A A 

One can understand these values as follows. By a rotation of —0(0) = — 7r/2 around the third 
axis the initial vector a(0) becomes (0, e, 7) T . Then by a rotation with angle —0(0) around 
the first axis it becomes (0, A, 0) T . Next, the rotation Ri(9(t)), followed by the rotation 
Rs((j)(t)) maps this fixed vector onto the time- dependent a(t). 
The Hamiltonian K(t) equals (see (63) of Part I) 

tsIx\ a 2^3 a CLid^ a\(X2 — 02^1 a 

ft\t) = — —bi — D 2 H 5 o 3 

Z A z z Z A 

= j|| [-4(t)S 1 + ai {t)a 2 {t)S 2 + a x {t)a z {t)S 3 ] + eS 3 . (16) 

The difference between this K(t) and the Hamiltonian H(t) as given by (jSJ) makes an extra 
rotation necessary. It involves the function a(t), given by (64) of Part I. It evaluates to 

a ( t ) = ai ^t(t) = 2eujkcn 2 (ut; k) 

U z 2 (ty {) ~ e 2 cn 2 (ujt;k)+u 2 sn 2 {ujt;k)dii 2 (u;t;k)' 1 ' 

The final result then becomes 

U(t) = e ^( i ) 5 3 e i(e(i)-e(0))5 le - i (7r/2)S 3e --i(/ t d S a( S ))( e 5i+753) > ( lg ) 

Note that this expression can be simplified to 

= e im)-7r/2)s 3e -i(0(t)-e(o))s 2e i\TX_ 

6jS* — \~ S 

with r = r(t) = — ds a(s) and X = . For further use note that 

Jo A 

e iXTX S je - iXrX = Sj + t sin(Ar)[X, Sj] + (cos(Ar) - 1)[X, [X, Sj}}. (20) 
2.4 General solution of the Bloch equations 

The general solution of the Bloch equations, given arbitrary initial conditions and with time- 
dependent Hamiltonian determined by (jHJ), can be derived using ( 1T91) . Note that (omitting 
time dependences and denoting 9q = 9(0)) 

2 2 

ecos(9-9 ) -7sin(fl-0 o ) = 6 ~ 7 cos(fl) - 2^ sin(0) 

A A 
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2 2 

jcos(9-9 ) + esm(9-9 ) = 2y cos(fl) + 6 ~ 7 sin(fl) 
7cos(#-# ) -esin(fl-# ) = -Asin(0) 

ecos(#-# ) +7sin(#-# ) = Acos(#). (21) 

These relations can be used to calculate 

U(t)(eS 1 + ^S 3 )U(t)* = -Xsm(9)S 3 + Xcos{9)S + , 
XU(t)leS 1 - -fS 3 )U(ty = [(e 2 - 7 2 ) cos(#) - 2e 7 sin(0) cos(Ar)] S+ 

+2 7 esin(Ar)S_ 

- [2e 7 cos(fl) cos(Ar) + (e 2 - 7 2 ) sin(0)] S 3 
U(t)S 2 U(t)* = - sin(Ar) cos(9)S 3 - cos(Ar),S_ - sin(0) sm(\r)S + , 

(22) 

with 

S + = sin(0)S'i + cos(0)S , 2 and = cos(0)S'i — sin(0)S , 2 . (23) 
Note further that 

(f> = — e + £ sin(0) tan(6>) 
= £cos(0) 
_ £sin(0) 

T " AcosW [M) 

Using these equations one can verify explicitly that the time-dependent operators ( 1221) satisfy 
indeed the von Neumann equation of motion. 

From (1221) one can obtain the three independent solutions a^ 2 \ of the Bloch 
equations. 

e / sin(0)cos(#) \ / sin(0)sin(#) \ / cos(</>) 

= — cos(0) cos(6 l ) — — cos(Ar) cos(</>) sin(6 l ) + — sin(Ar) — sin(0) 

A V -sin(0) J A V cos(#) / A V 

/ sin(0) \ / cos(0) \ / 

= — sin(#) sin(Ar) I cos(0) — cos(Ar) I — sin(0) — cos(#) sin(Ar) I 

V o ; v o / v 1 

/ sin(0)cos(#) \ e / sin(0)sin(#) \ £ / cos(</>) 

= — I cos(0) cos(6 l ) + — cos(Ar) I cos(0) sin(6 l ) — — sin(Ar) I — sin(0) 
A \ -sin(0) / A \ cos(0) / A V 

(25) 

The special solution ([9]) satisfies a = ea^ + 7 a ( - 3 - 1 . The general solution of the Bloch equations 
x = h x x is given by expression (67) of Part I and agrees with the results given above. 

2.5 The k = 1-limit 

We now focus on the case k = 1 because then the formulas simplify. The relevant expressions 
become 

sm{(f>) = — 

: V + c^tanh 2 ^) 



. N a 2 c<jtanh(a;t) 
cos(0) — 



Z \je 2 + c; 2 tanh 2 (^) 
a 3 wtanh 2 (u;t) - 7 

sm(0J 



cos(#) = 


z 
A 


= sech(. 


z = 


V 


^ a\ + a 2 . : 


A = 


V 


'z 2 + a§ = 


a = 




ai 

z 2 



U 2 + uj 2 tanh 2 (ujt) 



e 2 + i 2 



af + a 2 , = sech(u;t) \/ e 2 + u 2 tanh(u;i) , 



2ecu 



7 2 , 



e 2 + cj 2 tanh 2 (cjt) 

with 7 = (u; 2 — e 2 )/2u. The integral of a(t) can be done analytically. It yields 

r(t) = — [ dsa(s) = ^ let + arctan (— tanhftut)^ } . 
Jo e 2 + to 2 I V e J ) 

2.6 The Bloch equations at resonance 



(26) 



(27) 



The resonant condition is e = u. In that case, 7, as given by (TTOT) . vanishes. Hence, the 
normalisations z(t) and A, given by (|14|) and appearing in the expressions for the angles (f> 
and 6 (see (62) of Part I), become 



There follows 



z(t) = ujy cn 2 (cut; k)+ sn 2 (ujt;k)dn 2 (ujt;k), (28) 
A = uj. (29) 



sin(20) = = sn(2^t; k) (30) 

2 2 

cos(20) = fl2 ~ Ql = -cn(2a;t;fc) (31) 



2; 



sin(0) = = fcsn 2 (wt; fc), (32) 

A 



cos(fl) = -7 = \l cn 2 {ut- k) + sn 2 (ut; k) dn 2 (wt; k). (33) 
A " 

In particular is 9(0) = 0. The equations ( |22l) become 
U(t)S!U{ty = -sm(6)S 3 + cos{6)S + , 

U (t)S 3 U (£)* = sin(0) cos(c<jr)S' + — sin(u;r)S'_ + cos(u;r) cos(6 , )S , 3 
U{t)S 2 U{t)* = -sm{uT)cos{6)S 3 -cos{ujT)S_-sm{6)sm{ujT)S + . (34) 

They can be used to obtain the general solution of the Bloch equations ([6]) at resonance 

af n (t) = af n (O)cos(#)sin(0) 
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+ sin(0) sin(0) [a| en (0) cos(cur) - a| en (0) sin(cur)] 
— cos(0) [a| cn (0) cos(u;r) + a| en (0) sin(c<JT)] 

a| cn (t) = af n (O)cos(#)cos(0) 

+ sin(0) cos(0) [a| en (0) cos(^r) - af n (0) sin(^r)] 
+ sin(0) [a| en (0) cos(^r) + af n (0) sm{ur)) 

af n (t) = -af n (0)sin(#) 

+ [af n (0) cos(^r) - af en (0) sm{ur)) cos(9). 

The correction angle a, given by (IT7|) . simplifies to 

a (t) = = -ku (1 + cn(2o;t; k)) . 

z 2 

This expression can be integrated analytically. The result is 

dn(u)t; k) 



f 1 

r(t) = — dsa(s) = kt arctan 

Jo u 



k sn(a>t; k) cn(ut; k) 
With some effort, ([35]) can now be written as 



af a (t) = af n (0) cn(tut; k) 

+a| en (0) sn(u;t; k) cos(kout) 
+a| en (0) sn(u;t; k) sm(kujt), 
a f n (t) = af en (0) sn(W; k) dn(a>t; k) 

+af en (0) ( cn(a;t; k) dn(u;t; k) cos(kcot) — k sn(ut; k) sin(kujt)) 
+ a| en (0) ( cn(u;t; k) dn(a;£; k) sm{kojt) + k sn(W; k) cos(kujt)) , 
af n (t) = -af n (Q)ksii 2 (ujt;k) 

+a| en (0) (k cn(ujt; k) sn(u;t; k) cos(kut) + dn(u;t; k) sm(kut)) 
+ a| en (0) (k cn(ut; k) sn(a;t; k) sm.{kut) — dn(a;t; k) cos(ku>t)) . 



The three independent solutions are therefore 



cn(ujt; k) 
sn(u;t; k) d\i{uot\ k) 
—k sn 2 (a;t; k) 

sn(u;t; k) cos(kut) 
= ( cn(ut; k) dn{ujt; k) cos(kujt) — k sn(ut; k) sm(kuot) 
k cn(cot; k) sn(ut; k) cos(kut) + dn(a>i; k) sm(kcut) 
sn(ut; k) sin(fca;t) 
= ( cn(ujt; k) dn(u;t; k) sm{kujt) + k sn(u;i; k) cos(kut) 
k cn(cut; k) sn(a>t; k) sm{kuji) — dn(u;t; k) cos(kut) 

One clearly has a = ua^\ the special solution we started with. It is easy to verify that 
and a® are solutions of d = a x h. 
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3 SU(2) continued 

3.1 Including phase modulation 

A slightly different solution is obtained when the parameter e in is made time dependent 
in the following way 

e = e tanh(u;t). (40) 

In jl], the resulting equations are called the Bloch equations including phase modulation. 
Also in this case a special solution is known. The generalisation to Jacobi's elliptic functions, 
as given below, can be done in two different ways. 

3.2 Special solutions 

Assume that 

£(*) = £o™(ujt;k), 

e(t) = e sn(ujt;k). (41) 

or 

£(i) = £odn(ut;k), 

e(t) = e sn(ujt;k). (42) 

Both assumptions reduce to £(t) = £osech(u;t) and fj4*0|) in the limit k — 1. 
A solution of the equations ([6]) is given by 

ai(t) = eo cn(u;t; k), 

a>2(t) = w dn(u;£; k), 

a 3 (t) = -^ sn{ut;k). (43) 

or 

ai(t) = uj cn(cjt; k), 

a>2(t) = e dn(ujt;k), 

aM = -£ sn(ujt;k). (44) 

provided that £o = e o + uj2 ^ 2 '■> respectively £q = elk 2 + uj 2 , is satisfied. 

Only the first of the cases is treated below. The other case is completely analogous. 
Note that in the limit k — 1 the solution fl4*3l) reduces to the well-known solution 



ai(t) = eosech(u;£), 
a,2(t) = uj sech(ujt), 

a 3 (t) = -£ tanh(cL>t). (45) 



In the limit k = it reduces to a harmonic precession 

f(t) = £ cos(u;t), 
e(t) = e sin(o;i), 
ai(t) = e cos(u;t), 



a 2 (t) = uj, 

a 3 {t) = -£ sm(ut). (46) 
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3.3 Automorphisms 

Consider the special solution ( 1431) . The angles <fr(t) and 6(t) are determined by (62) of Part 
I. In particular, at t = is 

sin(0(O)) = |, cos(0(O)) = j, (47) 

and 6(0) = 0, with A = a/cq + lo 2 . This means that by a rotation i? 3 (— 0(0)) around the 
third axis the initial vector a(0) is rotated into the fixed vector A(0, 1, 0) T . Next, the rotation 
Ri(6(t)), followed by the rotation R 3 ((f)(t)) maps this fixed vector onto the time-dependent 
a(t). 

The Hamiltonian K(t) equals (see (63) of Part I) 

, . 0-2(23 „ aid^ 0\0 2 — 0-20,1 

Kit) = —bi-—S 2 + 6 3 

= j|| [-4(t)S 1 + ai (t)a 2 (t)S 2 + ai (t)a 3 {t)S 3 ] + eS 3 . (48) 

This is the same expression as (fT5l) . The difference between this K(t) and the Hamiltonian 
H(t) as given by ([SD makes an extra rotation necessary. It involves the function a(t), given 
by (64) of Part I. It evaluates to 

(A ai< ^ £(f\ - eo^o cnjut; k) dnjut; k) 

[) z 2 (t) X) e 2 cn 2 (ut;k)+u 2 dn 2 (ut;kY 1 J 

The final result then becomes 

U(t) = e^ ( * )53 e imSl e"^ (0)53 e iXT ^ t)x , (50) 

f>t ^ £j _j_ £j 

with r(t) = - / ds a(s) and X = — — . Note that (09} can be integrated analytically. 

Jo A 

The result is 

2u;AyU A — /i sn(u;t; k) 



with /i = \Jtl + k 2 
3.4 General solution including phase modulation 

It is now possible to obtain the general solution of the Bloch equations including phase modu- 
lation, with driving fields of the form f|4"Tl) . One calculates (again omitting time dependencies 
and denoting O = 0(0)) 

U(t)(e S 1 +uS 2 )U(t)* = \cos(6)S + - \sm(6)S 3 
\U{t){e S 1 - uS 2 )u{ty = [(e 2 - to 2 ) cos(#) + 2e o u;sm(0) sin(Ar)] S + 

+2eou cos(Ar)S'_ 

+ [2e uj cos(#) sin(Ar) - (e 2 - to 2 ) sin(0)] S 3 
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U(t)S 3 U(t)* = cos(Ar)sin(0),S + - sin(Ar)S_ + cos(Ar) cos(0)S 3 , (52) 

with 

S + = sin(<p)Si + cos(0)S , 2 and S- = cos(0)S'i — sin(0)5 , 2. (53) 

The three independent solutions of the Bloch equations including phase modulation are 
therefore 

sin(0) 008(6*) 
a«(t) = ^ I cos(0)cos(0) 



^ \ - Mi! 1 , 1 

w / sin(0)sin(0) \ ^ / cos(</>) 

+ — sin(Ar) I cos(0) sin(0) I + — cos(Ar) | — sin(</>) 



sin(0) 

i: 
»f 

cos(0) I A V 



(sin(0) cos(6 l ) 
cos(0) cos(0) 
- sin(0) 

(sin(0) sin(0) \ / cos(0) 

cos(0) sin(0) — - cos(Ar) I — sin(0) 
cos(0) / A V 

/ sin(0) ; 

a ( 3 )(t) = cos(Ar) cos(</>) sin(0) I - sin(Ar) I - sin(0) I. (54) 



sin(0) sin(0) ^ 




' cos(</>) \ 


cos(</>) sin(0) 


| " sin(Ar) | 


-sin(0) 


cos(0) j 




V / 



The special solution (l4"3"j) satisfies a = + ua^ 2 \ The general solution of the equations 

x = h x x has been given in Part I of the paper and can be rederived from the knowledge of 



a« a( 2 \ a( 3 ). 



3.5 The k = 1-limit 

In this limit the existence of the special solution requires that £o = A. The angle and the 
phase a become constants 

sin(0) = y and cos(0) = ^ (55) 
A A 

and a = — eo£o/A 2 . Hence one has Ar(t) = eot. The angle satisfies 

sin(0) = — tanh(u;t) and cos(0) = sech(co>t). (56) 

The three independent solutions are 

£q / e sech(u;t) \ ^ / -e tanh(u;t) \ ^ 

a^(t) = — I ujsech(ujt) J + — sin(e t) I — co>tanh(c<;t) J + — cos(e t) 

y Atanh(u;t) J \ Asech(cji) J 

(eo sech(u;t) \ / — e tanh(cjt) \ ^ 

u;sech(u;t) J — — sin(e t) I — u> tanh(cjt) J — — cos(e t) 

Atanh(cjt) J ' \ A sech^) J 
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j / — e tanh(u;£) \ ^ / u 

a (3) (t) = -cos(e t) -tutanh(cjt) --sin(e t) -e ]. (57) 
A \ Asech(c^) j X \ 

4 SU(1,1) 

The generators of SU(1,1) satisfy the commutation relations 

[S 1 ,S 2 ] = iS 3 , 
[S2,S 3 ] = -iSx, 

[S 3 ,S!] = -iS 2 . (58) 

Consider creation and annihilation operators and b satisfying the canonical commutation 
relations [b, b'] = I. Then the operators 

si = \m 2 +b% 

52 = \m*-b% 

5 3 = l -(tfb + btf) (59) 

satisfy (I58p . Hamiltonians which can be written as a linear combination of these generators 
are the generalised harmonic oscillators. 

4.1 Time-dependent frequency 

Consider the harmonic oscillator with arbitrary time- dependent frequency u(t) 

H = ^P 2 + l -muo 2 {t)Q 2 . (60) 

Here, Q is the position operator and P the momentum operator. They satisfy [Q, P] = i. 
This problem was studied by Lewis and Riesenfeld [2] . See also El Q31 IEI HSl El EH E21 EQJ - 
Introduce an annihilation operator defined by 

b ^2 Q + 'T2 P ' <61) 
with r = and ujq some constant frequency. Then the generators equal 

^ = ^ 2 -T p2 (62) 

5 2 = l(QP + PQ) (63) 

5 3 = ^Q 2 + jP 2 - (64) 
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Introduce the function j(t), modulating the frequency u , defined by u(t) = j{t)ujQ. The 
Hamiltonian becomes 

H(t) = -lu; (b-tf) 2 + ±u; 0l 2 (t)(b + tff 

= uo^ 2 {t) - l)Si + Ml 2 (t) + 1)S 3 . (65) 

Hence one has h = u; (7 2 (t) — 1,0, 7 2 (t) + 1) T . 
4.2 Special solution 

The time evolution equation reads (see the definition of the Lie bracket in Part I, equation 
(69)) 

a = h x a 

= (h 3 a 2 , h x a 3 - h 3 a x , /iia 2 ) T 

= w (7 2 (t) + l)(a 2 ,-a 1 ,0) T + ^ (7 2 (t)-l)(0,a 3 ,a 2 ) T . (66) 
Proposition 1 Let x(t) be a solution of the classical oscillator equation 

x + u 2 (t)x = 0. (67) 

Then a, defined by 

a = ^(i) 2 (l, 0, 1) T + ^ V(-1, 0, 1) T - lu xx(0, 1, 0) T , (68) 



is a solution of (E 

The proof is done by explicit calculation. One concludes that, to find a special solution 
of the von Neumann equation, it suffices to find a special solution of the classical equation 
( IBT|) . Note that the latter problem is equivalent with solving Riccati's equation 

<?-<? 2 =u; 2 (t). (69) 

The corresponding solution of (!67|) is 

x = Cexp ( - / dtg(t) J , (70) 



with integration constant C . In the case that u(t) is constant one finds 

g(t) = utan(ujt) (71) 
so that x(t) = C cos(u;£). With C = 1 and ujq = u> one obtains the special solution 

a = ^ 2 (- cos(2cut), sm(2ut), 1) T (72) 

If u(t) is of the form u{t) = ui (l + ecos(2At)) then the equation (J6"7|) is related to Mathieu's 

equation. The solution x = cn(u ,t; k), involving Jacobi's elliptic function with < k < 
1 / v2, is obtained when 



u{t) = lo ^2 dn 2 (co> t; k) - 1. (73) 
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4.3 Automorphisms 

Note that the special solution ( 1681) satisfies (a\a) = (using the metric with signature 
— ,—,+). Hence, we have to apply the exceptional case discussed in Part I, at the end of 
Section 5. 

Write a(t) into the form 

a(t) = i? 3 (0)Pi(x)i? 3 (-<KO))a(O). (74) 

From (77) of Part I follows 



CLi (x) 2 — UJqX 2 

a 3 (x) 2 + u 2 x 2 ' 



sin(0) 

a 2 2u xx 

COS(0) = — = — — g-j. (75) 

a 3 [x) z + unfix* 



From (83) of Part I follows 



a 3 (t) (if + uj-y 

x = ln M0j = l,, [(i)» + ^)]«- (76) 

Note that 0(0) simplifies if either x(0) = or i(0) = 0. 

From the general theory now follows that (see (76) of Part I) 

V(t) = e~ i(t> ^ S3 i e~ ix ® Sl e 1 ^ ^ 3 . (77) 

The corresponding Hamiltonian is 

K = x {cos(0)5i - sin(0)5 2 } + <j)S 3 . (78) 

Note that 

X — — = hi — and (p = 5 = n 3 — hi — . (79) 

a 3 a 3 a§ a 3 

Hence K can be written as 

K = (cos(0)x, - sin(0)x, 0) T 

1 T 

= — (hid 2 ,, — /i 1 a 1 a 2 , h 3 a\ — hiaia 3 ) . (80) 



One obtains h — k = aa, with 



a = hi — . (81) 



a 3 



The final result is 
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5 More general oscillators 



More general time-dependent harmonic oscillators have been considered in the literature 
[HI [151 l2Ql [221 EZ]- Even damped oscillators have been studied — see for instance 
[TT1 [TH [21]. Some of them can be treated by the present method. An example not yet 
considered in the literature is H = Y^j=i hj(t)Sj with 

hi(t) = acn(ut]k) (83) 
h 2 {t) = -asn(ut;k) (84) 
h 3 (t) = c + on dn(ut; k), (85) 

with constants u, a, c, and < k < 1. A special solution of the von Neumann equation is 
given by (see [25]) p s (t) = £^ =1 fl^t)^ witn 

a%(t) = aaa(ut;k) (86) 
02^) = — a sn(ut; k) (87) 
a 3 (t) = c. (88) 

Note that p s (t) is not a density operator. But this does not harm our method. 
Introduce unitary operators V(t) by (see (76) of Part I) 

V(t) = e - ^ 53 e~ ix ^ Sl e ix ^ Sl e^ ^ 3 (89) 

with angles <p(t) and x(t) satisfying (see (77) of Part I) 

sin(0) = cn(ut;k), cos(0) = sn(ujt; k) (90) 

C (X 

sinh(x) = , cosh(x) = ■ (91) 

Va — cr V« — c z 

Note that we assume that a 2 > c 2 . Also, note that 0(0) = ir/2. Because x(t) turns out to 
be independent of time t, (I8"9"j) simplifies to 

yty = e -Wt)-*(o))5 3> (92) 

Next calculate 

a(t) = = 1. (93) 

a z 

Hence, the time evolution is described by the unitary operators 

JJ(t) = V "(t)e~ ips(0) * = e-WityHtySa e -it(aSi+cS 3 ) _ 

It is now possible to calculate the time evolution of the generators in the Heisenberg picture. 
For simplicity take c = 0. Then one obtains 

Si(t) = u(tys 1 u(t) 

= cos(0(t) - 0(O))Si + sin(0(t) - 0(O))S 2 
14 



sin(0(£))5i — cos(0(t))S , 2 
cn(u;t; — sn(out; k)S 2 - 



(95) 



Similarly is 



u(tys 2 u(t) 



cosh(at) [cos(0(t) - (f>(0))S 2 - sin(0(t) - 0(O))£i] 
cosh(at) [sin(0(t))52 + cos(0(t))Si] — sinh(at)5 , 3 
cosh(at) [cn(a;t; fc)^ + sn(ut; k)Si] — sinh(at)<S , 3, 



sinh(at)S , 3 



(96) 



and 



U(t)*S 3 U(t) 

cosh(at)S 3 — sinh(at) [cos(0(t) — (f)(0)) S2 — sin(0(i) 
cosh(at)5 , 3 — sinh(at) [sin(0(i))S2 + cos(0(t))5i] 
cosh(at)5 , 3 — sinh(at) [cn(ut; k)S 2 + sn(a;t; k)S\] . 



(f>(0))Si) 



(97) 



Note that 5*3 is the energy of the unperturbed harmonic oscillator. Clearly, this quantity 
explodes for large times. Hence, the time- dependent harmonic oscillator described by (|85l) 
is at resonance. 

6 Discussion 

In Part I of this work a method was developed to solve time-dependent Hamiltonians with 
the assumption that they equal a time-dependent linear combination of generators of a finite 
Lie algebra. The method aims at finding all solutions, given one special solution. This 
Part II demonstrates how the method of Part I can be applied in two well-known cases, 
one corresponding with SU(2) symmetry, the other with SU(1,1). In this way, a number 
of known results are treated in a unified way. But in both cases the method is shown to 
produce new results as well. 

We did not try to reproduce the most general results found in the literature. We are 
confident that we could do so, at the expense of writing a more technical and less pedagogical 
paper. Of more interest is the application of our method to other Lie algebras. 

In [15], the Lie algebra that we used in Section 4 is extended to contain 6 elements |P 2 , 
\Q 2 1 \(PQ + QP)i P, Q, I- With this extension it becomes possible to calculate the time 
evolution in the Heisenberg picture of physically interesting quantities such as the position 
Q and the momentum P. One then can calculate the classical phase portrait by studying 
the orbit t — ► (ip\Q(t)\ifj) , (ifj\P(t)\ifj) for an arbitrary wavefunction ip. This will be done in 
a future work. 

Weigert [16] has considered the general case of SU(N) symmetry. Lopez and Suslov 
[28] used the HeisenbergWeyl group N(3) to describe a forced harmonic oscillator. Finally, 
Carinena et al [29], among others, are interested in developing a superposition principle for 
nonlinear equations by mapping the solutions onto the solutions of linear equations with 
time-dependent coefficients. 
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